Geometric frustration in the mixed layer pnictide oxides 
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We present results from a Monte Carlo investigation of a 
simple bilayer model with geometrically frustrated interac- 
tions similar to those found in the mixed layer pnictide ox- 
ides [Sr^MrizPn-iOi,Pn = As,Sb). Our model is composed 
of two inequivalent square lattices with nearest neighbor intra- 
and interlayer interactions. We find a ground state composed 
of two independent Neel ordered layers when the interlayer 
exchange is an order of magnitude weaker than the intralayer 
exchange, as suggested by experiment. We observe this result 
independent of the number of layers in our model. We find 
evidence for local orthogonal order between the layers, but it 
occurs in regions of parameter space that are not experimen- 
tally realized. We conclude that frustration caused by nearest 
neighbor interactions in the mixed layer pnictide oxides is not 
sufficient to explain the long-range orthogonal order that is 
observed experimentally, and that it is likely that other terms 
(e.g., local anisotropies) in the Hamiltonian are required to 
explain the magnetic behavior. 

I. INTRODUCTION 

Clean systems of interacting moments have been stud- 
ied extensively by analytic and numerical techniques. Al- 
though simplified models like Ising, Hciscnbcrg, and Hub- 
bard retain only the most fundamental interactions ob- 
served in real materials, they remain tractable to current 
theoretical techniques, and the study of their ordered 
phases in various regions of parameter space has con- 
tributed enormously to our understanding of magnetic 
phenomena and the physics of correlated systems. ||] 
However, real materials are never clean. There is often 
frustration due to competing interactions and disorder in 
the interaction strengths. 

Competing interactions that cause magnetic frustra- 
tion can have many origins, lattice geometry, magnetic 
and non-magnetic impurities. In three dimensions, he- 
lical magnetic order has been observed when geomet- 
ric frustration is accompanied by anisotropy. |^,|| Spin 
glasses phases are observed when frustration is accom- 
panied by random disorder. ||-^| It has also been sug- 
gested that some non-collinear spin ordered structures 



belong to a new chiral universality class. [[7| The systems 
we study are essentially two dimensional and contain no 
anisotropic terms or disorder. Frustration is caused by 
the lattice geometry. Our primary focus is the orthogonal 
magnetic structure observed in the mixed layer pnictide 
oxides. H 

The pnictide oxides of type A2Mn 3 Pri202{A = 
Ba, Sr; Pn = As, Sb) are layered antiferromagnets that 
contain two distinct square planes of manganese atoms 
arranged in a lattice of space group symmetry IA/mmm. 
In one layer, manganese is bonded to oxygen in a pla- 
nar CuOi arrangement, MnO^ . In a second layer, it 
is bonded to a pnictogen in a tetrahedral structure, 
MnPn 2 ~ , where pnictogen atoms project alternately 
above and below the plane defined by the manganese 
atoms. From here on we denote the two layers as Mn(l) 
for MnOf~ and Mn(2) for MnPnj~ . The Mn atoms 
from the two planes are arranged so that a site in the 
Mn(l) layer sits directly above and below the center of a 
square plaquette of Mn atoms in the Mn(2) layer. The 
manganese carry a spin S — 5/2. Frustration can enter 
through nearest neighbor interlayer coupling. A more 
detailed investigation of these systems has been reported 
elsewhere. M 

II. MODEL 

In all the pnictide oxides except Sr2Mn^As202 there 
is long range order in the planes that eventually gives 
rise to weak ZD order. In the compound SrvMnsAszO^, 
there is only short range order in the Mn(l) planes. 
The ability of ordered planes to drive c-axis order has 
been investigate before in the case of layered antiferro- 
magnets. ph| |l2]] In Sr2Mn 3 Sb2C>2, magnetic order in 
the Mn(l) layers is established along the a-axis of the 
magnetic unit cell, while in the Mn(2) layers the mag- 
netization is along the c-axis. Hence there is an or- 
thogonal alignment between neighboring layers. Such 
an ordered state is not without precedent. Jl3|,[l4| How- 
ever, the different temperatures at which the layers order 
(Tm„2 ~ 300A" and TmuI ~ 65A) and the symmetry of 
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the frustrated interlayer interactions which leads to can- 
cellation suggest a system of two independent Neel or- 
dered layers. Previous work has shown that this is not 
always the thermal or quantum fluctuations (in 

frustrated systems) can lift the degeneracy of the system 
to select a single state. |l5| , [l^1 

To study the effect of frustration on the ground state 
magnetic order of the pnictide oxides, we developed a 
simple model of classical Heisenberg spins with nearest 
neighbor intra- and interlayer interactions. The basic 
structural unit is a set of two layers, one each of type 
Mn(l) and Mn2 (see Fig. 0). The Mn(2) layer has a 
lattice constant a = 1 and contains n 2 sites. The Mn(l) 
layer is larger by a factor \[2 and is rotated by 7r/4 with 
respect to the lattice directions of the other layer. The 
Mn(l) layer contains n 2 /2 + n + 1 spins. Note that the 
interlayer coordination is not the same for spins on the 
two layers. A spin on the Mn(l) plane is coupled to four 
spins on the Mn(2) plane, and each spin in the Mn(2) 
plane is coupled to only two spins on the Mn(l) plane. 

The Hamiltonian for our bilayer model is written as 



H = Jl y sf } • s {l \ +j 2 y • s { % 

i,<5i z,(52 



are over nearest 
one has \S\ = 



The constants J\,J%, and J± represent the Mn(l) and 
Mn(2) intralayer couplings and the interlayer coupling, 
respectively. The summations of 5 P 
neighbors to site i. For classical spins 
(S 2 X + S 2 + S 2 Z Y' 2 = 1. The relatively large spin-5/2 of 
the Mn atoms in the pnictide oxides makes this a reason- 
able approximation. 

We studied the equilibrium physics of our model by 
a single spin flip Monte Carlo algorithm. We have ad- 
dressed concerns about proper sampling of phase space 
by performing simulations with random and ordered ini- 
tial configurations. We have also considered the effects 
of the boundary on our finite simulations by employing 
a few different boundary conditions: open, periodic, and 
periodic with an effective field on the Mn(l) edge sites. 
In all cases considered, we found no qualitative differ- 
ence in our results due to the initial configuration or the 
conditions imposed at the boundary. 

To determine the relative orientation between neigh- 
boring spins, either within the same layer or in different 
layers, we measured a collinear 
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and a perpendicular 



cr 13 = 
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spin-spin correlation function. Here summations are per- 
formed over all nearest neighbors S of site i and then over 
all sites in the lattice; z is the coordination number and 
N a is the number of sites in layer a. Intralayer correla- 
tions are denoted by a = (3 and interlayer correlations are 
represented by a ^ (3. We stress that C\\ and C± measure 
local correlations. For classical Heisenberg spins, these 
correlations take on the simple forms Cii = (cos 2 6) and 
C±_ = (sin 2 #). In the high temperature, paramagnetic, 
limit, the values C\\ — 1/3 and C± = 2/3 are obtained. 
We also measured the magnetization and staggered mag- 
netization of each layer. 
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FIG. 1. A 2D projection of the two distinct layers of the 
pnictide oxide Sr2MnsSb202- Sites in the Mn(l) layer are 
represented by dark circles while sites in the Mn(2) layer 
are represented by light circles. The intralayer couplings are 
shown as Ji and J2, and the interlayer interaction is indicated 
byJx. 



III. RESULTS AND DISCUSSION 

From the experimental data, representative couplings 
would set the Mn(2) intralayer exchange to be stronger 
than the Mn(l) intralayer exchange, with the interlayer 
interaction weaker by at least an order of magnitude. 
Therefore, experimentally motivated couplings in our 
model were set to J2 — 2.0, J\ = 1.0 and J± = 0.1. For 
the results presented herein, systems with 1600 spins per 
Mn(l) layer and 841 spins per Mn(l) were equilibrated 
for 15, 000 to 25, 000 sweeps followed by 15, 000 measure- 
ment sweeps with 10 — 25 sweeps between measurements. 
One Monte Carlo sweep denotes an update of all spins 
on the lattice. 

Our results for a bilayer model indicate that, as a func- 
tion of temperature, the moments within each layer be- 
gan to order when the temperature dropped below the 
respective energy scale, e.g., T = J\ for Mn(l) moments 
and T = J2 for Mn(2) moments. However, the eventual 
ground state was a system of two Neel ordered layers 
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with an arbitrary orientation between the magnetization 
directions. In a simulation with four layers (i.e., a lay- 
ered sequence Mnl — Mn2 — Mnl — Mn2) and periodic 
boundary conditions along the c-axis, we observed the 
same qualitative behavior as a function of temperature 
(see Fig. ||). We emphasize that the interlayer spin-spin 
correlation function, C,, ' , remained at the paramagnetic 

limit down to low temperatures, T ~ J±. 
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FIG. 2. Temperature dependence of the local intra- and in- 
terlayer spin-spin correlations in the four layer model with pe- 
riodic boundary conditions and Ji = 1.0, J2 = 2.0, J± = 0.1. 
A parallel alignment is favored for intralayer spins when the 
temperature drops below the respective intralayer coupling; 
however, the interlayer correlations remain at the high tem- 
perature limit of 1/3 even for T ~ J±. 

We also studied the effect of the strength of frustra- 
tion on the magnetic ground state. To do this, we fixed 
the temperature and swept in values of Jj_. For a bi- 
layer model with J± < 0.25, we observed two Neel or- 
dered layers with a paramagnetic interlayer orientation, 
i.e., Cm' ~ 1/3 or independent layers. As Jj_ was in- 
creased, the two Neel ordered layers moved towards a 
collinear state. We observed this behavior independent 
of the initial configuration (see Fig. ||). 

An explanation of these results can be understood by 
considering the physics of two simpler models, a zigzag 
lattice and a diagonal lattice, (refer to Ref. §). In the 
zigzag lattice the frustration is along one crystalline di- 
rection; this geometry is similar to the Mn(2) to Mn(l) 
interlayer interaction. In the limit of weak Jj_, a uni- 
formly canted state is established in the zigzag lattice 
with orthogonal order obtained in the limit Jj_ — > 0. In 
the diagonal lattice, one has the equivalent of the J\ — J 2 
model. This interlayer geometry is analogous to the cou- 
pling of Mn(l) sites to Mn(2) sites. In the limit of weak 
Ji , the diagonal lattice is composed of two Neel ordered 



layers with a collinear alignment. Hence, at weak frustra- 
tion, it is the competition between these two tendencies 
in the experimental model that leads to a paramagnetic 
orientation with large fluctuations. In the case of strong 
frustration, J± > 0.5, the zigzag and diagonal geometries 
both select a collinear interlayer arrangement. 

An orthogonal state for our bilayer model can be found 
but this phase occurs at couplings strengths that are not 
supported by experiments (refer to Ref. ||). In the case 
where J2 = 0, the resultant model is a network of inter- 
secting zigzag chains. By setting J± — 1 and sweeping in 
Ji a transition to a uniformly canted state was observed 
in the limit of large J\. For a model with J2 = J± = 1.0, 
we observed a Neel ordered Mn(2) layer and a paramag- 
netic Mn(l) layer, but with Mn(l) spins orthogonal to the 
local Mn(2) environment, for J x < 0.25. At J% > 0.25 the 
system tended toward a state with collinear alignment. 

We conclude that frustration caused by nearest neigh- 
bor interactions, both intra- and interlayer, in the mixed 
layer pnictide oxides is not sufficient to explain the long 
range orthogonal order that is observed experimentally. 
However, recent work on mixed metal pnictide oxides 
does indicate that the two distinct layers do not behave 
independently. C^] In these systems, it is likely that other 
terms in the Hamiltonian, e.g., ion anisotropics arising 
from spin-orbit effects, are required to explain the mag- 
netic behavior. 
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FIG. 3. Interlayer spin-spin correlations as a function of 
J± with Ji = 1.0 and J2 — 2.0. The results are for a bilayer 
model with the initial configuration of the layers being either 
random or Neel ordered. 
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